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We study the transverse single spin asymmetry (SSA), An, for the single inclusive 
process £p^ — > h + X, in a perturbative QCD factorization scheme with inclusion 
of spin and transverse momentum dependent (TMD) distributions. By adopting the 
relevant TMDs (Sivers and Collins functions) as extracted from semi-inclusive deep 
inelastic scattering (SIDIS) and e + e~ data, predictions for these SSAs are given. A 
measurement of An for this process could then provide a direct test of the validity of 
the TMD factorization. 

1 Introduction 

In this contribution [1] we discuss a preliminary study of transverse SSAs in inclusive hadron 
production in lepton-proton collisions as a tool to test the TMD factorization hypothesis. 

Azimuthal and transverse single spin asymmetries, with their behaviour and size, defi- 
nitely represent a challenge for the QCD factorization theorems and, at the same time, have 
opened a new way to learn on the spin structure of hadrons. A large amount of data has been 
analyzed by many experimental collaborations in various inclusive processes (for reviews see, 
e.g., Ref. [2]): the left-right asymmetries in single polarized proton- proton collisions (from 
the early E704 to the latest STAR and BRAHMS data), the azimuthal asymmetries in 
single polarized SIDIS (HERMES and COMPASS) and in hadron-pair production in e+e~ 
annihilation (Belle). 

Two theoretical approaches have been proposed to describe such phenomena [2]: i) a 
generalization of the pQCD factorization scheme with inclusion of a new class of spin and 
TMD distribution and fragmentation functions; ii) a collincar QCD factorization theorem 
in terms of higher- twist parton correlation functions. 

The first approach was initially adopted [3] to explain the early SSAs observed in pp — > 
tt+X and later applied, with success, to predict [4] An at much larger energies. An extended 
and systematic study can be found in Ref. [5] . Presently, the TMD approach is believed to 
hold for SSAs characterized by the presence of two scales, a large (i.e. Q 2 ) and a small one 
(pt — Aqcd), like for the production of small px hadrons in SIDIS processes or of small pr 
lepton-pairs in Drell-Yan processes [6, 7]. For a detailed and complete classification of all 
TMDs and their role in azimuthal and SSAs in SIDIS see Ref. [8]. 
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Based on these results the first extractions of the Sivers, Collins and transversity func- 
tions have been performed [9, 10, 11, 12, 13, 14, 15]. Notice that all these analyses of SIDIS 
data are carried out in the 7* — p cm. frame, according to the following TMD factorization 
formula: 

da^' h+x ~ £]/«,/,,(*, *j.;Q a ) ® do*-*** ® D h/g (z,px]Q 2 ) , (1) 
1 

where fc^ and p± are, respectively, the transverse momentum of the quark in the proton 
and of the final hadron with respect to the fragmenting quark. 

The alternative formalism, the higher-twist approach, has been proved to hold for SSAs 
where a single hard scale is relevant, like the inclusive production of large pt hadrons in 
hadron-hadron collisions [16]. A corresponding and rich phenomenology has been developed 
in Rcfs. [17, 18, 19]. 

In the last years these two approaches have been shown to be somewhat related and 
equivalent in the kinematic regime where they both apply [20] . 

However, a definite proof of the validity of the TMD factorization for SSAs in inclusive 
particle production in hadron-hadron collisions, with a single large scale, is still lacking. 
In this context we mention that a sort of modified TMD factorization approach has been 
discussed in the study of dijet production at large px in pp collisions by including the proper 
gauge-links structure in the elementary processes [21, 22, 23]. 

What we aim at discussing here is then an experimental test of the TMD factorization 
by considering the SSA in the £p — > h + X process, where a single large pt final particle 
is detected [24] . This is the analogue of the left-right asymmetry observed in pp — > h + X . 
Differently from SIDIS, we will now perform the analysis in the lepton-proton cm. frame 
without the detection of the final lepton, but still requiring a large Q 2 regime (see below). 

Similar studies, although with different motivations, were presented in Refs. [25, 26, 27]. 

2 TMD approach to £p — > h + X 
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We consider the process yfil — > h + X, with p moving 
along the positive Z-eods in the p — £ cm. frame, the final 
hadron produced in the X — Z plane and | direction along 
+Y (Fig. 1). We assume the same TMD factorization 
scheme as for the process p^p — > h + X and compute 
Am = (da 1 * — da^)/{da^ + da^). The numerator of this 
SSA reads 

da 1 - da 1 ~ ^{a^/^t cos^ <g> da <g> D h/q Figure 1: Kinematics. 
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dAa <g> A N D h/q i cos 0c + hf^ /p <g> dAa <g> A N D h/q7 cos(</> c - 2tf> q )\ , 

(2) 



where the q£ —> q£ elementary cross sections are given by 



2 s 2 +u 2 2 su 

da oc e„ s — ctAcr oc — e„ — , (3) 
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with s,t,u the Mandelstam invariants and <j>c = ^ + <Pq', with 4> q {(f) q i) the azimuthal angle 
of the initial (fragmenting) quark in the cm. frame and <f>^ the hadron azimuthal angle 



in the outgoing-quark helicity frame. Notice that the result shown in Eq. (2) can be also 
obtained as a particular case of the (A, Sa) + {B, Sb) —^C + X process [5]. In equation (2) 
we recognize the contributions from the Sivers function [28, 29] (first line) 

f g /pt (x, fej_) - fq/pi (x, k±) = A N f g / p i (x, k±) S-fpxfeJ, (4) 

where the mixed product in our configuration gives the cos <p q dependence in Eq. (2), and 
from the Collins function [30] (second line) 

D h/g r(z,pi_) - D h/ql {z,p±) = A N D h/g r(z,pi_)s q ■ {p q x p ± ) , (5) 

coupled with hi (the TMD transversity distribution) and hj~ T ; in this case the azimuthal 
dependences in Eq. (2) arise from the phases entering the TMD distributions, the elementary 
scattering and the Collins function. 

Before presenting our results we want to comment on some aspects peculiar for this 
process. First of all, as in SIDIS, we have a single partonic subprocess and only the t- 
channel contribution (much simpler than the pp — > h + X case). However, without the 
detection of the final lepton one is not able to reconstruct the lepton plane and then to 
access separately the Sivers and the Collins effect (as usually done for SIDIS). On the other 
hand in the backward region (w.r.t. the proton direction) |w| becomes smaller and so does the 
spin transfer dAa/da (see Eq. (3)) implying a strong dynamical suppression of the Collins 
effect. At the same time, this does not affect the Sivers contribution since, contrary to what 
happens in pp — ■> h + X, no w-channel is active . [Remember also that the t variable strongly 
depends on <p q (the azimuthal dependence of the Sivers effect, see Eq. (2))]. 

It is also well known that the ip — > h + X process is dominated by quasi-real photon- 
exchange with low pt hadrons in the final state. However, for the validity of the adopted 
perturbative QCD factorization scheme the elementary scattering must be governed by a 
large momentum transfer, let us say Q 2 > 1 GeV 2 . This is trivially guaranteed in the usual 
collincar configuration by requiring a moderate-large pt for the final hadron (j>t > 1 GeV). 
We have checked that, by including k± effects in the kinematics, the safe region (Q 2 > 1 
GeV 2 ) corresponds to a minimum pr value of the order of 1.5 GeV. More precisely, for 
such value one has to restrict to the backward region (xf = 2p^j^fs < ), while, at larger 
values, like pr = 2.5 GeV, we can explore the full xf range, remaining in the large Q 2 
regime. These are the values considered in our estimates. 

A word of caution has still to be added. In our calculation we do not include any hard 
gluon emission which could give a large pr hadron even in the low Q 2 regime. This would 
lead to a two-jet event and might be experimentally excluded by requiring the absence of 
any hadron activity in the opposite hemisphere w.r.t. the detected hadron. 

3 Results 

We show our estimates for two different kinematic setups related to the ongoing HERMES 
(PLab = 27.5 GeV) and COMPASS (pLab = 160 GeV) experiments and focus on pion 
production. By numerical calculation we have checked that the (maximized) contribution to 
An coming from hj~ T (see Eq. (2)) is totally negligible, as well as the analogous contribution 
to the unpolarized cross section involving the Boer-Mulders function. In our estimates we 
will adopt the present information on the relevant TMDs (Sivers, Collins and transversity 
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Figure 2: A N {pH -> n + X) for HERMES kinematics (p Lah = 27.5 GeV) at fixed p T values 
as a function of xf- Left and central panels: Sivers effect (the statistical uncertainty bands 
for charged pions at pr = 1.5 GeV are also shown); right panel: Collins effect. 
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Figure 3: Same as in Fig. 2 but for COMPASS kinematics (p Lab = 160 GeV). 



functions) as extracted in the latest analyses of SIDIS and e + e data [14, 15]. We will show 
separately the contribution to Am from the Sivers or the Collins effect alone. 

In Figure 2 (left and central panels) we present our estimates for the Sivers contribution 
to Am at HERMES kinematics for two values of pr- As discussed above, for pt = 1.5 
GeV only the backward region can be considered. For charged pion production at pt = 1.5 
GeV we also show the statistical uncertainty bands from the fits. The largest Am values 
obtained correspond to the x region (in the polarized proton) where the Sivers functions, for 
up and down quarks, reach their maxima. Notice that while for positive xp, the minimum 
of x (in the polarized proton) is given, roughly, by Xf, for negative xf the minimum of x 
is controlled by the ratio px/y/s, with z always bigger than \xf\- It is interesting to note 
the sizable Am for ir~ production (larger than the corresponding Aut in SIDIS) due to the 
dominance of the down-quark contribution with a small contamination from the up quark. 

For pt = 1.5 GeV the Collins effect (not shown), involving hi, is negligible, while it 
reaches at most 1-2% at the largest \xf\ values for pr = 2.5 GeV (Fig. 2, right panel). 

In Figure 3 we show the analogous results for COMPASS kinematics. Again at pr = 1.5 
GeV only the Sivers effect gives a sizable contribution (left panel), while the Collins effect 
(not shown) is compatible with zero. At pt = 2.5 GeV the Sivers effect (central panel) 
dominates only in the backward region, while in the forward region the Collins effect (right 
panel) becomes sizable. 

Notice, however, that for xf > 0.3 we are probing the Sivers and the transversity func- 
tions in a region (x > 0.3) where they are not constrained by present SIDIS data. 



For larger energy values, e.g. y/s = 100 GeV, like those reachable at the proposed EIC 
experiments and at the same px values as above, we obtain negligible contributions to the 
SSAs from both the Sivers and the Collins effect. This is due to the low x region explored 
and the corresponding sea-quark dominance. 

Conclusions 

We have presented a preliminary phcnomenological study of SSAs in Ip — > h + X within a 
TMD factorization scheme. This process presents interesting aspects somewhat in between 
the SIDIS case (for which the factorization is generally accepted) and the pp — > h + X case 
(where it is assumed) . In this sense it could represent a clear test of the TMD factorization 
hypothesis. By adopting the parameterizations so far extracted from SIDIS data, sizable 
SSAs (~ 5 — 10%) can be obtained, mainly due to the Sivers effect and, to a lesser extent, 
to the Collins effect. 
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